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^ t Abstract 

Large discrepancies have been observed between measured Electromagnetic 
Dissociation(ED) cross sections and the predictions of the semiclassical Weizacker- 
Wilhams- Fermi (WWF) method. In this paper, the validity of the semiclassical 
5_j ■ approximation is examined. The total cross section for electromagnetic excita- 

tion of a nuclear target by a spinless projectile is calculated in first Born approx- 
imation, neglecting recoil. The final result is expressed in terms of correlation 
functions and convoluted densities in configuration space. The result agrees 
with the WWF approximation to leading order(unretarded electric dipole ap- 
proximation) , but the method allows an analytic evaluation of the cutoff, which 
is determined by the details of the electric dipole transition charge density. Us- 
ing the Goldhaber- Teller model of that density, and uniform charge densities 
for both projectile and target, the cutoff is determined for the total cross sec- 
tion in the nonrelativistic limit, and found to be smaller than values currently 
used for ED calculations. 

In addition, cross sections are calculated using a phenomenological mo- 
mentum space cutoff designed to model final state interactions. For moderate 
projectile energies, the calculated ED cross section is found to be smaller than 
the semiclassical result, in qualitative agreement with experiment. 



1 Introduction 



The availability of relativistic heavy ion beams has opened a new avenue for the 
study of electromagnetic excitations of nuclei. Cross sections are enhanced both by 
the charge of the projectile ions, and by the relativistic contraction of the projectile's 
electric field into a sharp pulse of radiation at high energies. Experiments range from 
single and double nucleon-removal reactions Jl]], to the study of "halo" nuclei using ra- 
dioactive beams0, to the possibility of multi-phonon excitations of collective nuclear 
states||. Aside from their intrinsic interest, electromagnetic excitation processes in 
peripheral collisions will also be important at RHIC|4]]. 

There is, therefore, a tremendous incentive to develop an understanding of the 
physics involved in these processes. Such an understanding has two facets. On the 
one hand, one must be able to calculate the detailed structure of the target nucleus 
in order to calculate its response to a particular probe, such as the electromagnetic 
interaction. Fortunately, this has been one of the central topics of nuclear physics 
for many years, and an extensive literature exists on the subject ||. On the other 
hand, it is also necessary to understand the process by which the projectile excites 
the target. With few exceptions to date||, this aspect of the problem has been dealt 
with by using the venerable Weizacker- Williams-Fermi (WWF) method of virtual 
quanta @, and its generalization to arbitrary multipoles||. This approach is based 
on the observation, due to Fermi, that the electromagnetic fields of a point charge, 
when boosted to high energy, are transverse to the direction of the charge's motion. 
The supposition is that one can calculate cross sections by replacing the projectile 
by an equivalent pulse of electromagnetic radiation. The cross-section for a given 
reaction is then 

cr ww = J dojn{uj) cr 7 (co>), (1) 

where cr 7 (u;) is the cross section for the same reaction induced by real photons, and 
n{uj) is the number of photons of energy u in the pulse of equivalent radiation. For 
dipole transitions, n(u) can be determined by calculating the intensity of the projectile 
fields as a function of frequency, integrated over impact parameters^]. Assuming the 
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projectile's trajectory is a straight line, the result is 

n El {u) = -^(xKoWK^x) - ^rX 2 [Kl{x) - K*(x)}) 

2Z„ 2 a 1.123, (3\ 
~ — log )-^r) as 7^00, 2 

where i^o(i) is a Bessel function of imaginary argument, Z p is the projectile charge, f3 

7/3 



its speed, x = U3hm f ; n , and 6 m j n is an impact parameter cutoff required to get a finite 



result. 

In some respects, the Weizacker-Williams-Fermi(WWF) method is very well suited 
to the problem of electromagnetic dissociation by heavy ions. Cross sections may be 
calculated either using a model, or, as is often done, by direct appeal to measured 
photodisintegration cross sections. In the latter case one obtains a cross section, at 
least at high energy 0], that is independent of any model of the target's structure. 
Furthermore, the equivalent photon number, n(uj), can be obtained by a relatively 
straightforward classical calculation. Nonetheless, the WWF method is not com- 
pletely satisfactory. Theoretically, there is no systematic procedure for evaluating 
corrections to the semiclassical result, and consequently no way of gauging the relia- 
bility of the approach. In addition, there is the troublesome problem of choosing the 
minimum impact parameter. Since 6 m , n is not fixed by the WWF procedure itself, 
a number of choices have been proposed [11]], based on phenomenological considera- 



tions. On the experimental side, large discrepancies from the naive WWF predictions 
have been observed in target fragmentation experiments [12|], leading some authors to 



question the validity of first order perturbation theory for heavy projectiles JHJ. 

The aim of this paper is to go beyond the semiclassical WWF method and calcu- 
late the quantum mechanical cross section for electromagnetic excitation of nuclei in 
heavy ion collisions. In the next section, we derive an expression, originally used in 



high energy physics |L4] , for the unpolarized cross-section in first Born approximation, 
neglecting recoil effects for both the projectile and target nuclei. We demonstrate that 
the cross section does approach the WWF approximation in the limit of large projec- 
tile energy, and in an appendix we derive an expression for the cutoff parameter, 6 m j n , 



for dipole transitions, using a simple model for the transition matrix elements. In the 
following section, we use the same simple model to compare our results for the full 
cross section to those of the WWF approximation as a function of projectile energy, 
transition multipolarity, and transition frequency. Finally, using a phenomenological 
cutoff designed to model final state interaction effects, we compare the results of this 
simple model with measured single-neutron removal cross sections. 

2 Quantum Excitation Cross Section 

We begin this section with a discussion of the kinematics of the nuclear excitation 
process. The projectile nucleus, with mass Mi and momentum Pf 1 = (Ei, Pi), scatters 
from the target by exchanging a virtual photon of momentum g M . In the process, 
depicted in figure 1, the target of mass M 2 is excited from its ground state to an excited 
state of mass M2+00. For high Z p projectiles, elastic scattering of the projectile(oc ZT) 
will dominate over inelastic processes(cx Z p ), so that in the following we shall assume 
that the projectile remains unexcited in the final state. This is in accord with the 
semiclassical picture, where the projectile remains in its ground state, following a 
straight line trajectory throughout the collision. Kinematically, the elastic scattering 
of the projectile translates into a condition on the four momentum transfer q, 

q 2 = 2P-q - q 2 = 2EMo-P-q), (3) 

— » — » 

where (3 = Pi/Ei is the projectile velocity. For nuclear transitions, the momentum 
and energy transfers are on the order of tens of MeV, while the projectile and target 
masses are on the order of tens of GeV. As a result, q 2 /2Ei is negligible, and we 
immediately obtain the minimum momentum transfer 

2 <?o 

\q\rrUn = qo/P " > = ( 4 ) 

where 7 = —j= ^ . In exactly the same fashion, we derive a relation between the 
energy transfer and the target excitation energy, 

q = uj + O{uj 2 /M 2 ). (5) 
4 



Even at this stage of the calculation, we see a similarity with the semiclassical calcu- 
lation at high energy. For large projectile energies, the minimum momentum transfer 
goes to zero, so that it is reasonable to expect that the cross section will be dominated 
by the pole in the photon propagator at q 2 = 0, where the interaction cross section 
will be well approximated by the cross section for real photons. 

Having treated the kinematics, we now turn to the calculation of the cross section. 
The method we use is a straightforward evaluation of the Feynman diagram of figure 
1, borrowing heavily from the theory of relativistic electron scattering 15| in order 



to relate the target and projectile form factors to the appropriate nuclear transition 
matrix elements. For a + projectile (the generalization to nonzero spin is trivial), 
the spin averaged cross section for an arbitrary transition multipole is given by 

(^z pa ) 2 f d 4 q d 3 P f , ^ 4jt 4, 



aw =- ((W-S?W a J (2^2E f (2ny {27T)6{Pf - Pl + q) 

xF;(q 2 )PtP?WT u (q 2 , q -K t ), (6) 

where the matrix element of the projectile current is given by Z p F p (q 2 ) (P i/Lt + 
and the spin-averaged target structure function is given by 

d 3 K f ,_ 1 



(2E Kf )(2n)^ ' y 1 '(24 + 1) 

Y,{KiMJM\ K f, m ')( K fMJ»{V>)\Ki,™) (7) 

mm' 

= W 1 (q 2 , q .K i )(g, v - q -^) 

q A 

+W 2 (q 2 ,q-K l )(K t ,- q -^)(K w - q -^), (8) 

where m(m') denote the azimuthal angular momentum quantum numbers of the ini- 
tial(final) state, (Ex f ,Kf) is the final state four momentum, and W\(q 2 ,q ■ Kj) and 
W2{q 2 , q ■ Kj) are Lorentz invariant target structure functions. For what follows, it is 
critical to realize that while the structure functions Wi(q 2 , q ■ Ki) and W 2 (q 2 , q ■ Kj) 



are different for different transition multipoles, the tensor form of W^ v is determined 
by current conservation and parity, and thus independent of both the multipolarity 
of the transition and spin of the target. Using the elastic scattering condition for the 
projectile, and keeping only the leading terms in inverse powers of the projectile and 
target masses, we get 

x MfW^q 2 , q ■ Ki) + (Pi ■ K^W^q 2 , q ■ K,)] . (9) 

At this point, it is useful to reexpress the two invariant structure functions in terms 
nuclear matrix elements. In particular, we shall choose the two rotationally invariant 
combinations (Jo(O)Jo(O)) = (p(O)p(O)) and (J(0) ■ J(0)), where the brackets are in- 
troduced as a convenient shorthand for the spin sums and matrix elements of equation 

7. Hence, in the target rest frame, 

r cfiK 



^ -<f(J(0)J(0)) + q i(p(0)p(0)) 

x w { } 



;; (2g 2 + g 2 )(p(0)pW)-g^(J(0)J(0)) (n) 

For analytical purposes, it is convenient to have expressions in configuration space. 
Replacing the momentum delta functions by an integral representation, and using 
translational invariance, we obtain 



x I Sz c m -f(J(S)A0)) + (12) 

J 
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x r dh e m W + gg)MfM<>j> - fiMMl ( i3) 

J q 
The projectile form factor is related to the projectile rest frame charge density via 

F(q 2 ) = f d 3 xe^'-*p p (x), (14) 



where q ' = q — ^t-/3 ■ q/3 is the three momentum transfer in the projectile rest frame. 

Since the matrix elements appearing here are rotational scalars, we may replace 
the complex exponential by its angular average with impunity. The result for the 
cross section then becomes, after performing the integrals over the three momentum 
transfer, 



(q°yj {2kY2E K} 
[ d 3 z[(p c (z)p c (0) - J c (z)J c (0))(^I 1 (q z,P) + I 2 (q z,P)) 



+(p c (z)p c (0))(^h(q O z,f3) + I 3 (q z,f3))} , (15) 



where 



/•oo 
i qdq- 



jojqy) 
l - ?)< 



/•oo 
l qdq 



j-[(Si((l/P - l)y) - Si((l/P + l)y)(cosy + y siny) 
+(Ci((l//3 + l)y) + Ci((l/(3 - l)y))(siny - y cosy) + 2/3 7 2 sin(y//3)], 

(16) 

Jo (92/) 



g 2 (l - g 2 ) 



= I[i oosy(5t((l/)9 - l)y) - + l)y)) 

y 
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+ | siny(«((l//3 + l)y) + Ci((l/(3 - l)y)) + /3 sin(y//3) - yCi(y/P)], 



(17) 



^ 3 (y,/3) = 3 



■oc 



—3oi.lv) 



W 2 - y 2 mj (y/P) + t- cos(y//3) + ^ C%//3), 



(18) 



with ^(y) and Ci(|/) the sine and cosine integral functions, respectively, and 



At this juncture, the expression bears little resemblance to the compact semi- 
classical result in equations 1 and 2. There is no simple factorization of the integrand 
into a flux factor multiplying the photo-cross section. The integrand is not even a 
function of the same variable, q /^(3, as the expressions in equation 1. Asymptotically, 
the cross section of equation 14 falls off like a power for large q°, in contrast to the 
exponential decrease dictated by the Bessel functions of the semi classical result. In 
light of this, it seems quite unlikely that the semiclassical expression will yield a good 
approximation for the cross section over a large energy range. 

On the other hand, the correspondence principle tells us that the two expressions 
must agree in the classical limit. For the problem at hand, the classical regime 
occurs for large projectile energies, where the straight line trajectory of the projectile 
represents a good approximation to the classical trajectory. In that limit (7 — > 00), 
the momentum integrals become 





(19) 



h(y,P) 



jjo(y) + Oil) 



(20) 





(21) 



h(y,P) 



(22) 



and the cross section is given by 
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CT\VW 



{Z p af f dq° 



d 3 K f 



(2ir)5{E Kf -M 2 -q c 



x J d 3 z[(p c (z)p c (0) - J c (z) ■ J c (0))Uq°z) (1 - log( 7 2 )) + 0(1)]. (23) 

Using methods identical to those just described, the cross section for real photons 
can be written as, 

ira f d 3 Kf 



q° J (2tt) 3 2E K/ 



{2n)5{E Kf -M 2 -q l 



x J d 3 zj (q°z){p(z)p(0) - J(z)J(0)). 



From which it follows that 

2(Z p a) 2 



(Tww 



dq° 



(q 



0\2 



log( 



27 



q°R/ lE 2 



(24) 



+0(l/ 7 2 



(25) 



where all the terms of order (I/7) have been absorbed into the cutoff parameter Ri. 

The first term in this expression is just the high energy limit of the semi-classical 
cross section, with no cutoff parameter. The cutoff parameter, Re, represents the 
leading order corrections due to Coulomb and longitudinally polarized virtual pho- 
tons, as well as off shell corrections to the photonuclear cross section. In the limit of 
large 7, the precise value of log(q°Ri) will be negligible compared to log(j), and the 
semiclassical result is recovered. 

In the low frequency limit, the off-shell corrections to the photonuclear cross sec- 
tion vanish, and Re is determined solely by the projectile and target transition densi- 
ties. In Appendix B, we exploit this fact to calculate R\, using the Goldhaber- Teller 
model [jl6|| for the transition densities. Accurate determination of this parameter, 
which governs the total excitation cross section, will be relevant for estimating the 
large background due to electromagnetic processes at RHIC, where the semiclassical 
limit should be realized. 
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3 Results 



In this section, we compare results both with the semiclassical calculations using the 
Weizacker- Williams method, and, ultimately, with data from single neutron removal 
experiments. To do this, we are required to assume a model for the transition matrix 
elements appearing in eq. 14. For simplicity's sake, we again choose the Goldhaber- 
Teller(GT) |TB| model to describe the transition densities, and assume that the pro- 
jectile and target are described by uniform density spheres. A brief sketch of the 
GT model, as applied to dipole excitations, may be found in the appendix. For gen- 
eral multipoles, a more detailed description of the model, including electromagnetic 
transition matrix elements, may be found in ref. 15. 

In order to effect a comparison with the semiclassical calculation for arbitrary 7, 
it is useful to decompose the total cross section into multipoles. Since the magnetic 
multipole contributions to the cross section are small, we shall restrict our attention 
to electric multipoles larger than zero. The decomposition may be accomplished by 
inserting a factor 

1 = — J d 3 xd 3 x'S 3 (z-x + x'), (26) 
into eq. 14, and by using the Bessel function identity 

J (q\x-x'\) = ^-J2Uqx)Uqx')Y em (x)Y; m (x'), (27) 

Ira 

to do the angular integrals. The result is 

%w = ^§^E / W-Wi) f qdqF*( q >) 

rj (Ffa) - F*j(g)) + (i?(g)-Fj(g)) + 5M + 1 Fib) I 
L7 2 (qi-q 2 ) 2 {ql~q 2 )q 2 q A i 2 (ql~q 2 )q 2 ^ 

(28) 

where, for t > 0, F l p (q) = £C e j 2 (qR t ), and Fj(q) = ^±ll^F e p (q), with R t the radius 
of the target nucleus and Ct a constant chosen such that the photo-cross section 
satisfies the energy weighted sum rule for multipole £, and F p (q') = 3ji(q'R p )/q'R p is 
the elastic form factor of the projectile, with R p the projectile radius. 
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For comparison, we also calculate the semiclassical cross section, using a prescrip- 
tion from reference 11 for the minimum impact parameter, 



brain 1.34 



Al /3 + a] /3 - o.75(a; 1/3 + a; 1/3 )\ , (29) 



appropriate for single and double nucleon removal experiments. The ratio of the 
quantum and semiclassical cross sections for 12 C and 197 Au projectiles incident on a 
197 Au target are shown, as a function of projectile energy, in figures 2 and 3, respec- 
tively, for both El and E2 transitions. We assume that the E1(E2) transition is to 
a sharp, isovector(isoscalar) giant resonance state at 13.8(11.0) MeV, and that the 
energy weighted sum rule is saturated. We further assume that the giant resonance 
state decays exclusively via one neutron emission. For high projectile energies, we 
find that the quantum El cross section is enhanced by about 10 per cent relative to 
the classical result. This enhancement agrees with the results of the appendix, where 
it was shown that the quantum mechanical cutoff is smaller than that of eq. 29, 
resulting in a larger cross section. More surprising is the enhancement of the E2 cross 
section, which remains a factor of two larger than the classical result at 7=100. Also 
noteworthy is the huge enhancement of both the quantum El and E2 cross sections 
at small projectile energies. For 7 >5, nearly all of the difference in the El cross 
section may be reabsorbed into a redefinition of b m i n as described in the appendix, so 
we conclude that the WWF method is a good approximation for the total El cross 
section for energies of 5 GeV/nucleon and higher, provided that the minimum impact 
parameter is properly chosen. For lower energies, much of the difference can be elim- 
inated by altering b min , but discrepancies as large as 20 per cent persist at very low 
projectile energies. 

In figures 4 and 5, the ratio of the quantum and semiclassical cross sections is 
shown as a function of energy for 7=2, assuming that the momentum dependence of 
the transition matrix elements does not vary with transition energy. As u goes to 
zero, the El cross section approaches the semiclassical result. This is expected, as in 
this limit the argument of the logarithm is large in both the semiclassical expression 
and the high energy limit of the quantum calculation. As a result, the semiclassical 
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piece of the quantum cross-section dominates, and the two limits are the same. For 
large u, both the El and E2 ratios are enhanced, and the degree of enhancement 
is quite sensitive to the behavior of the transition form factors. For small to, the 
quantum E2 cross section is roughly four times the classical result, indicating that 
the large 7 and small u limits are not equivalent for higher multipoles. 

In order to make a meaningful comparison with experiment, we must account for 
the final state strong interactions between the target and projectile. The choice of 
b min in eq. 29 is designed to do exactly this for the semi-classical problem |ll]]. Briefly, 



b min is chosen such that the mean number of nucleon-nucleon collisions, calculated 
using the Glauber approach , at impact parameter b = b m i n is exactly one. For 
b > b m i n , the probability that a second nucleon will be knocked out of the target 
by the strong interaction drops rapidly to zero. For b < b m i n , the mean number 
of nucleon-nucleon collisions rises very rapidly, and the probability of at least one 
additional nucleon getting out of the target goes rapidly to one. Hence, we obtain 
the usual semiclassical picture, where the equivalent photon number is calculated by 
integrating the field intensity from all trajectories with b > b m i n . 

The same idea can be applied to the quantum picture by artificially including a 
maximum value, q m ax ~ 1/^mm, for the transverse momentum transfer in the collision. 
The effect of this cutoff is to change the upper limit of the q integration in eq. 28 
to \J ql/ f3 2 + \jb 2 min instead of 00. In figures 2-5, the ratio of the new quantum 
and semiclassical cross sections is shown as a function of the projectile energy and 
transition frequency for the same projectile target combinations used previously. The 
same general trends hold; the quantum cross section is enhanced for small projectile 
energies, and approaches the semiclassical limit as 7 gets large. (This latter fact holds 
because the cross section is dominated by the photon pole as 7 —>■ 00, so that the 
precise value of the high q cutoff becomes irrelevant.) Not unexpectedly, the size of 
the quantum cross section is smaller than in the case of no cutoff, with the result 
that the small 7 enhancement of the cross section is less pronounced, and, for large 
7, the semiclassical limit is approached from below, rather than above as before. 
Similar conclusions may be drawn regarding the transition frequency dependence of 
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the cross section. As before, the small u and large 7 limits of the El cross section 
approach the classical result, while the two limits differ markedly for the quadrupole 
cross section. Perhaps the most noteworthy feature of the uo dependence is that the 
strong dependence of the cross section ratio on the transition form factors has been 
largely eliminated. 

Of particular phenomenological interest is the 20-25 per cent suppression of the 
quantum El cross section relative to the semiclassical in the region 7 = 2 — 3. This 
suppression is at precisely the right location and magnitude to explain the discrep- 
ancy between the WWF approximation and recent single neutron removal data for 
238 U on 197 Au|12[. In table 1, we compare the results of the our calculation and the 
WWF calculation with experimental data from references 2 and 12. For the heavi- 
est projectiles, the agreement between the quantum theory and experiment is much 
improved over the WWF results. For Fe and Ar, the quantum theory does about as 
well as the semiclassical, and for the lightest two projectiles, the semiclassical theory 
does better. 

This level of agreement with the quantum theory is in fact quite heartening. To 
begin with, the "experimental" numbers listed in table 1 are not raw data, but rather 
the difference between the raw data and an estimate of the strong interaction contri- 
bution to the single neutron removal cross section. As noted in reference 11, estimates 
of the strong interaction contribution to the cross section are model dependent, so 
that an additional systematic uncertainty of at least 20 mb should be added to all 
the results listed in table 1. Clearly, for the heavy projectiles, this additional un- 
certainty is not important, as the extracted electromagnetic cross sections are quite 
large by comparison. For the lightest projectiles, however, the additional uncertainty 
is comparable to the extracted ED cross section, and is quite likely responsible for 
the observed discrepancies. 
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4 Conclusions 



Using a very simple model of nuclear structure, we have calculated the electromagnetic 
excitation cross section for nuclear collisions in first Born approximation, neglecting 
final state interactions, and have found significant differences from results of the 
semiclassical WWF method. For El transitions, we find that the impact parameter 
cutoff required for the semiclassical calculation to agree with the high energy limit of 
the quantum cross section is significantly smaller than the phenomenological cutoffs 
used to analyze experiments. At energies less than a few GeV/nucleon, the cross 
section is enhanced over the semiclassical result by as much as 20 per cent when the 
smaller cutoff is used. For E2 transitions, we find the cross section is significantly 
enhanced even at RHIC energies, and, unlike the El case, that the limit of small 
transition energy and large projectile energy are not the same. 

When final state interactions are included via a phenomenological cutoff, we find 
that the El cross section is greatly enhanced over the analogous WWF calculation 
for low projectile energies, while the cross section is suppressed at higher energies. 
The E2 cross section is suppressed by as much as factor of three for all but the 
lowest projectile energies. The low energy enhancement is relevant for electromagnetic 
dissociation studies of the low energy, pygmy resonances in neutron rich nuclei such 
as n Li, and the suppression at higher energies resolves the conflict between single 
neutron removal experiments and the semiclassical theory for all but the lightest 
projectiles, where the data is very sensitive to systematic errors in the separation of 
the nuclear and electromagnetic contributions to the cross section. 

While the agreement with data is extremely satisfactory given the simplicity of 
the Goldhaber- Teller model, a number of interesting questions remain. To begin, 
the calculation should be redone using a realistic model for the projectile and target 
densities in order to improve the quantitative description of the data. In addition, 
the effect of additional photon exchanges should be studied, both to understand the 
effect of the repulsive coulomb potential on the scattering process, and to study the 
question of multi-phonon excitations in the target nucleus. Finally, the process where 
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both the target and projectile are excited, which should be non-negligible in target 
fragmentation experiments with light projectiles, should be calculated. 
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Appendix A 

In this appendix we sketch the cross section derivation for target excitation in a 
more familiar form, which highlights the nonrelativistic nature of the nuclear physics. 
It also allows us to point out our approximations (explicit or implicit) and to empha- 
size the role of gauge invariance, which puts the results into a form where Siegert's 
theorem can be applied immediately. We assume the following: 

1. First Born approximation in the fine-structure constant, a; this is required for 
tractability. 

2. No target recoil, which presupposes that momentum transfers are very small 
compared to the target mass, m t ] this greatly simplifies the kinematics. 

3. The maximum momentum transfer can be replaced by infinity, for ease of per- 
forming integrals; in practice, small momentum transfers dominate. 

4. A nonrelativistic description of the target nucleus, and of the internal structure 
of the projectile nucleus, is sufficient; photonuclear physics is almost entirely 
nonrelativistic, and most of what we know is based on this (successful) descrip- 
tion. 

5. Both target and projectile are spinless and only the former is excited; neither 
restriction is essential, but will simplify the derivation. 

6. All purely hadronic contributions to the excitation amplitude can be ignored. 
Elastic Coulomb scattering is infinite and will also be ignored. 

The current of the elastically scattered projectile is conserved and is given by 
F p (q 2 )(Pf + PiY / \J 'AEfEi, where the projectile charge distribution is present through 
its Fourier transform, the projectile form factor, F p (q 2 ). The latter is a function 
of the (squared) four-momentum transfer, g M = P i — Pj. The distinction between 
this quantity and the three-momentum transfer is a relativistic correction (and a 
recoil correction in the projectile rest frame). We consequently replace F p (q 2 ) by 
F p (c?) = Jd 3 xp p (x) e^ s . 
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The target current is denoted by J^(q) = (p(q)), J(q)) and is conserved. That is, 
g M J M (g) = 0. Rewriting the kinematical factors in the projectile current, (P/ + p) M , as 
2Pf — g M , we can use current conservation and drop the q^ factor. Thus, the current- 

1 — * — * — * — * 

current matrix element is given by (Ei/Ef)* (pno(q) — P • Jno(q)), where (3 = Pi/Ei 
is the projectile velocity and J^ (q) =< N\J^(q)\0 >. The energy difference go of the 
states labelled N and is denoted by un- 

We wish to calculate the total cross section, which means that we must integrate 
over all of phase space. Ignoring recoil (which means that Pj is independent of the 
scattering angle, 9), we have for fixed energy transfer: q 2 = (Pf — P) 2 = P| + Pf — 
2PfPiCos6, or dq 2 = 2PfPisin6d8 = dQ(2PfPi/2ir), since the azimuthal dependence 
is trivial. All of the phase-space integrals (but one) can now be performed and we 
obtain: 

°™ = ~I,i.(F* l '" ,,v ' 8 ' J "' i ' (A1) 

This deceptively simple form has many aspects of complexity. The integral should be 
dominated by the pole in the photon propagator and the small values of q^ min = uo 2 N / f3 2 
and q^in = —uf^ / f3 2, y 2 . This further implies that electric dipole processes should 
dominate, because they are largest for small q 2 . This multipole should consequently 
not be very sensitive to F p , unlike higher ones which are certain to be. In addition, we 
can make use of the conserved current to transform to Coulomb gauge. This eliminates 
the (redundant) longitudinal component of the current, and expresses the result in 
terms of transverse-current matrix elements (which determine photoabsorption) and 
purely Coulombic excitation. We also impose the no-recoil approximation in the form 
Pj • q/Ei = q Q — (3 ■ q = q 2 /2Ei ~ 0, so that we can replace (3 ■ q by uo N (= q ) where 
needed. We also impose current conservation: q • Jno(q) — Pno(q)- 

Expanding the squares of matrix elements and using the two spin-summed rela- 
tionships 

E PnMJ%M = 4 E \pm(q)\ 2 , (A2) 

spins " spins 
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E JnM4*M= E [^\pnM?^ p + \JnM? - ^\r^iv: 



3u 2 



spins 



spins 



X 



(5 af3 - q a f) 



leads to the (spin-summed) result: 

\pN0 ~ P ■ JnO? = \pNo(q)\ 2 



7$ , .2 \2 



00 



NJ 



2 V" <f , 



(A3) 



oo 



\Jm(q)\ -^-\PNo(q)\ 



(A4) 



Note that the factor which multiplies (the first) |p| 2 cancels the photon propagator 
and substitutes the Coulomb one. Moreover, the quantity in square brackets is the 
transverse current (i.e., the full current minus the longitudinal part). This leads to 
the Coulomb gauge result: 



0"WW 



\pm(q)\' 



+ 



UP 



V - oo 2 -V 



N) 



\JNM\ 2 -^\Pmm 



{Ah) 



Using (3 2 = — 4j + 1 this can be rearranged into a form commensurate with eq. (14), 
after we perform the ^-integrals. We use 

J^ (g) = J d 3 xe^ s < N\j»(x)\0 > (A6) 

and the fact that there is no overall dependence on q after spin sums are performed, 



which leads to a slightly different form of eq. (14): 



4ttZV 



<7ww — 



P 2 



E ~T f d * x f d 3x> \pno(x)p* no (x) 



— + h(y) 



Jno(x) ■ Jno(x) - Pm{x)p* m {x' 



T 
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(A7) 

One projectile form factor has been folded into each charge and current matrix el- 
ement. That is, Pno(%) is actually the convolution of a projectile (elastic) charge 
density with a target (transition) density. We have also used y = oj^z = u^\x — x'\. 
This compact result in configuration space is finite, requires no cutoffs, and is rela- 
tively simple to work with. 

The corresponding expression for the cross section for the absorption of a photon 
with energy, u, can be developed in the same form: 

a 7 (uj) = 2n 2 a> ^ — — — — f d 3 x [ d 3 x' j (y) \j m (x) ■ J* m {x* ) - p m (x)p* m {tf) . 

(A8) 

Various sum rules can be constructed from this by integration: 

a n = J duo LU n a~ f (uj) . (A9) 

In Appendix B the leading-order contributions are worked out in detail, including 
an analytic evaluation of the cutoff, it^(for I = 1, in the context of a simple nuclear 
model. 



20 



Appendix B 

Complete expressions including formulae for the cutoffs Ri are difficult to develop. 
They will also differ from multipole to multipole. Nevertheless, because most of 
the photodisintegration cross section at low energies is unretarded electric dipole in 
nature, it is useful to investigate how the cutoff arises and estimate its size according 
to models of the electric-dipole transition density. 

If one expands to lowest order in z (viz., (z)°) the factors multiplying the current 
matrix elements [Jjvo(x) • Jno(^')} m ec f (A7), only unretarded electric-dipole tran- 
sitions can result. Using Siegert's theorem ||15|| to express the current in terms of the 
dipole operator, D, and the excitation energy, u N , 

J d 3 xJ N0 {x) = iu N D m , (Bl) 

one finds that the current terms contribute to the total disintegration cross section 

AnZ 2 a 2 

cww = — — T N0 , (B2) 



an amount 

4 J o } = lAvo| 2 (ln( 7 ) - y) • (B3) 

Note that these terms neither require nor generate a cutoff in the logarithm. 

The corresponding constant factors multiplying p N0 (x)p* m (x') give no net contri- 
bution to this order, since 

Jd 3 xp N0 (x) = Q m = , (B4) 

because we have agreed not to include elastic scattering (N ^ 0). The first-order 
terms in z 2 generate elastic scattering plus electric-dipole transitions: 

TnI = 1 [d 3 x [ d 3 x'p m (x)p* m (x')z 2 \n(z) 



| |Srf ( ln (-) + i ln(7)+7E _^_a, (B5) 
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where 7^ is Euler's constant, z = x — x', and z 2 = (x 2 + x' 2 ) — 2x ■ x! has been used. 
Combining the charge and current contributions one obtains 

1 



T, 



NO 



d x I d x p N0 (x)p m (x ) z \n(z) 



+ ^\D N o\ 



In 



^7 



- IE 
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~6~ 



This can be rearranged into the conventional Weizacker- Williams form 

2 



T 



NO 



\D 



NO 



\ \uj N Ri 



-IE 



if we define: 



ln(iV2) = -— - f d3x f d 3x 'P m ^)p*No(x') z 2 \n(z) 
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2\D m \ 2 



(B6) 



(B7) 



:B8) 



Given any (electric-dipole) transition density, the p-wave part of z 2 \n(z) can easily 
be projected out and the double integral performed. The cutoff,i?i, comes from the 
J 3 -term in eq. (A7), which itself comes solely from Coulomb scattering (i.e., the first 
term in eq. (A5)). 

If the integration variables in eq. (B8) are changed to z and x', we find 

d 3 x / d 3 x' pno(x)p* no (x') z 2 \n(z) 



d 3 zF N (z)z 2 \n(z) 



(B9) 



where 



F N (z) = J Sx p m {z + x)p* m (x) , (BIO) 
which is a convolution of transition densities. Elastic densities of this type have 



appeared in atomic calculations |Tj, [18]. This form is particularly convenient for 
deriving asymptotic forms, as well. The drawback is the required construction of the 
second-order (convoluted) transition density. We note that the logarithm in eq. (B8) 
is the only one contained in Is(y), but it contributes to all Coulomb multipoles. 
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A complete model of the electric-dipole transition density is needed in order to 
construct F^{z). One such model is the Goldhaber- Teller model |TBJ, which we sketch 



below. The transition charge operator pno(x) can be represented by the surface- 
peaked function 

p m (x) = (Nm\p{x)\0) = —Xno Kn ■ Vpo(x) , (Bll) 

where (0|p(x)|0) = po(x) is the ground-state charge density, from which it follows 
that 

J d 3 xxp m {x) = \ m e* m = D N0 . (B12) 

We have defined / d 3 x p (x) = 1 and used the spherical projection operators e* m and 
the Wigner-Eckart fll| theorem to determine the dependence of the transition matrix 



element on the magnetic quantum number (m) of the final (dipole) state. For this 
model then, |Z?tvo| 2 — l^ivol 2 an d the last term in eq. (B8) can be rewritten as 

d 3 xp {x) J d 3 x'p (x')(ln(z) + -) 



= — + - J d 3 xpo(x) J d 3 x'po(x') [(x + x') 2 \n(x + x') — (x — x') 2 In \x — x'\]/xx' 
where the elastic counterpart of F^(z) is the Zemach density [jT7f|: 



p (2) (z) = J d 3 x' Po (\z + x'\)p (\x'\) . (B14) 

The latter density plays a critical role in the nuclear-size modification of the hyperfine 
splittingfl7| and Lamb shift |H| in atoms. Note that if po is properly normalized, then 



P(2) is also. The final form in eq. (B13) results from performing the average over the 
angle between x and x' contained in z. This collective model then leads to the very 
simple result: 

ln(i2i/2) = -1 + / d 3 zp {2) (z) \n(z) . (B15) 



23 



To go further requires another model assumption. For simplicity we choose a 
uniform (liquid-drop) charge density: po( r ) — ^^j^ , which leads to 



The integral in eq. (B15) has the value ln(2i? t ) — |, which leads finally to 



(B16) 



ln(i? 1 /2) = ln(2^) - ? - , (B17) 



or §l = 0.695, where R t = 1.2 A? fm. 

This derivation ignores one important piece of physics. The projectile, as well as 
the target, has a finite size. Because in momentum space the form factors of initial 
and final projectile and target are all multiplied together, the projectile charge density 
will enter as a convolution with the (transition) density of the target. For example, 
folding the target density in the previous model (with radius Rt) with a projectile 
whose density is p' Q (with radius R p ) leads to a transition density 

p(x) = J d 3 x'p' (\x-x*\)? m -V , p (x , ) = t m -Vp {2) (x) , (B18) 

where p( 2 ) is obtained by folding p and p' Q together. Thus, we can use eq. (B15) 
if we fold two densities (p(2)) together. This density, p( 4 )(x), can be shown to be 
independent of the order in which the folding is performed. Consequently, one can 
form p( 4 ) by first folding two densities po together, then separately folding together 
two densities p' , and finally folding these two (folded) densities together. 

The resulting forms are rather complicated, and we simply state the result: 

(]n(z)) = J d 3 zp (4) (z) ln(z) = ln(2(R t + R p )) 



^ 8 (R P + 2R t )(R t + 2R p )(2Rl + llR t R p + 2Rp 
+ ^- R ^ 2100RJRJ 



Rt — Rp 



Rt + Rp 
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^ (4^-81^) ^ / R p \ _ (4g-8Lg) ^ f R, 



p 1050i? t 6 \R t + Rp) 1 1050i?6 \Rt + R P 



R A p R\ 107R 2 t 107g 16661 



525Rf 525R* 700R 2 p 700R 2 12600 



There are two simple limits: 



and 



Um Rp ^o [(ln(z)> - ln(2(^ + R p ))] = -~ 

Zzm^ [<ln(z)> - ln(2( J R t + i^))] = -H ln(2) - ^ . (B20) 

The terms in eq. (B19) which supplement \n(2(R t + R p )) (which we denote by s(^-)) 
vary monotonically from -0.75 to -1.114 as R p is varied from to R t . Thus we can 
rewrite eq. (B15) as 

\n(R 1 /2) = \n(2(R t + R p ))-l + s , (B21) 

or 

Ri = S (R t + Rp) , (B22) 

i 

where R Pjt = 1.2 A p t , and S = Ae~ l+S varies from 0.70 to 0.48 for the range of 
variation of R p above. Previously, a value of 1.0 was recommended 0. 

Although the model used here can be criticized on the basis that it is too simple to 
be realistic, our derivation nevertheless represents an exact solution to the problem 
for a simple collective model of electric-dipole transitions (Goldhaber- Teller model 
with a uniform charge distribution) in the unretarded dipole limit. The troublesome 
cutoff is determined by details of the electric-dipole transition charge density. 

Finally, we can express our results in a particularly simple and elegant form which 
emphasizes well-known photonuclear sum rules. We find 

a ^ = —^ 2 -{ ln {^R 1 )- lE -Y) ' (B23) 

4 

where <r_i is defined in eq. (A9) and has the experimental value, 0.22(2) Af mb for 
medium to heavy nuclei, and uj is the mean photon energy, which can be taken to be 
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the giant resonance energy, 79 A t 3 MeV |[20|| . The factors in front of the bracket are 

4 

(numerically) consistent with Z^Af / (3 2 fib. The argument of the logarithmn is £/?7, 
where £ varies from 6.0 to 4.3 as R p varies from to R t . 
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Figure Captions 
Figure 1 Electromagnetic excitation process. 

Figure 2 Ratio of the quantum to classical excitation cross sections for 12 C on 
197 Au as a function of projectile energy. 

Figure 3 Ratio of the quantum to semiclassical excitation cross sections for 
197 Au on 197 Au as a function of projectile energy. 

Figure 4 Ratio of the quantum to classical excitation cross sections for 12 C on 
197 Au as a function of transition energy. 

Figure 5 Ratio of the quantum to classical excitation cross sections for 12 C on 
197 Au as a function of transition energy. 



29 



a class( m ty a quantum^ mh ^ <J EXPT (mb) 

Proj 7 El E2 El E2 



12 
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C 3.23 42 9 36 4 75 ± 14 



Ne 3.23 111 22 95 9 153 ± 18 

40 Ar 2.91 315 62 262 25 348 ± 34 

56 Fe 2.81 614 120 506 47 601 ± 54 

139 La 2.34 2,190 462 1,738 171 1,970 ± 130 

238 U 2.0 4,337 1,045 3,388 365 3,160 ± 230 



Table 1 Predictions of Weissacker- Williams- Fermi and quantum theory for single neu- 
tron removal from 197 Au using the Goldhaber- Teller model for the transition densities. 
Data from references 2 and 12. 
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